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The origin of charge inhomogeneity in YBa2Cu30i / is investigated using a new experimental 
method designed to determine the nuclear quadrupole resonance (NQR) asymmetry parameter r\ 
for very wide NQR lines at different positions on the line. The method is based on the measurement 
of the echo intensity as a function of the angle between the radio frequency field Hi and the principal 
axis of the electric field gradient. Static charge inhomogeneity deduced from 77 > are found in this 
compound, but only in conjunction with oxygen deficiency. This limits considerably the possible 
forms of charge inhomogeneity in bulk YBa2Cu3 !/ . 



The discussion on the mechanism for high tempera- 
ture superconductivity (HTSC) is focused these days on 
the presence or absence of charge and spin inhomogene- 
ity in the CuC>2 planes. Such inhomogeneity can lead 
to a one dimensional boundary, possibly in the form of 
stripes, between "hole-rich" and "hole-poor" regions and, 
allegedly, to superconductivity. pj Indeed, it is now estab- 
lished by both surface and bulk techniques that most of 
the underdoped cuprates phase separate into "hole-rich" 
and "hole-poor" regions. For example, there is a con- 
sensus derived from surface sensitive STM experiments 
on underdoped Bi2Sr2CaCu20s+(5 that the planes are 
inhomogeneous. A different example is La2- a; Sr a ;Cui04 
[LSCO] where the evidence for phase separation is de- 
rived from experiments sensitive to both magnetic fluc- 
tuation, such as muon spin relaxation (/iSR) Q, and 
charge fluctuations such as nuclear quadrupole reso- 
nance (NQR).0] Even in very underdoped YBa2Cu30 y 
(YBCOy) phase separation is observed; neutron scat- 
tering from phonons related to charge inhomogeneity is 
found in doping levels up to YBCO6.35 and A*SR de- 
tects a spin glass phase for a similar doping range. 
However, the origin of this phase separation is still not 
clear; is it coming from competing phases in the Cu02 
plane, or does it simply stem from the quenched disorder 
introduced by doping? 

A possible way to address this question is to apply 
NQR to YBCOy since in this compound one can dis- 
tinguish between different in-plane coppers [Cu(2)] reso- 
nance lines, and associate each line with a local environ- 
ment. However, a standard NQR based determination 
of charge homogeneity is not available to date due to 
the complicated and very wide spectrum of underdoped 
YBCO (which will be reviewed below). The purpose of 
this work is to overcome the NQR problems and to shed 
light on the evolution of charge homogeneity in the bulk 
of YBCO close to optimal doping. For this purpose we 
developed a new experimental technique based on the 
measurement of the Cu(2) NQR echo intensity as a func- 
tion of the angle between the radio frequency (RF) field, 
Hi, and the principal axis of the electric field gradient 
(EFG) at the copper site. We call this angle-dependent 



NQR technique ADNQR. 




FIG. 1: (Color online) A toy model calculation for 1 77 1 in the 
presence of charge stripes in the CUO2 plane. See text for 
details. 



Standard NQR is based on the fact that nuclei with 
spin I > 1/2 can be viewed as positively charged oval 
objects. As a result, their energy inside a solid depends 
on their orientation in the electrostatic potential V(r) 
generated by the other nuclear and electronic charges. 
When the nuclear poles are close to positive charges their 
energy is high, and when the poles are close to negative 
charges the energy is low. The energy difference between 
different orientations is determined by the EFG tensor 
Vij = q x .q x at the position of the Cu nuclei. The di- 
rections can be chosen so that Vij is diagonal. These 
directions are known as the principal axis of the EFG. 
Due to Laplace equation (V xx + V yy + V zz = 0) the NQR 
Hamiltonian is determined by only two parameters, v q 
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and 77, and is given by 

n q = ^[3i 2 z -i 2 + v(i*-i 2 y )] (i) 

where v q is a frequency scale proportional to the EFG 
V zz and 



is a dimensionless number. It is customary to choose the 
directions so that \V XX \ < \V yy \ < \V ZZ \, and therefore 
< j] < 1. For the spin 3/2 Cu nuclei this Hamiltonian 
has only one resonance frequency given by 

/ = + (3) 

In YBCO7 and YBC06 the directions are known ex- 
perimentally: z is the c direction and x and y are di- 
rections in the CuC>2 plane 0, Q- 111 YBCO6.5 these 
directions are determined to be the same on theoretical 
grounds and are believed not to vary for all other 
doping 9]. Assuming the directions are doping indepen- 
dent, r\ is a measure of the local charge homogeneity of 
the CuC>2 planes. When these planes are homogeneous 
with local xy rotation symmetry 77 = 0. In contrast, 
when the planes are inhomogeneous due to phase separa- 
tion as in the case of static stripes, for example, then on 
the boundary between hole poor and hole rich stripes the 
xy rotation symmetry will be lost and we expect r\ ^ 0. 
This situation is demonstrated in Fig. ^ obtained using 
a toy model. The parameters of this model are chosen 
for clarity only and have no implications on the physical 
situation in YBCO. A plane of a square lattice with total 
charge Q is sandwiched between two similar planes with 
total charge —Q/2 each. The electric field produced by 
each ion is screened with a screening length of two lat- 
tice sites. The charge in the central plane is distributed 
in the form of stripes, 10 atoms wide, and rj is numerically 
calculated. In the figure only the central plane with its 
stripes, and |?7| along one line of ions crossing different 
stripes, are shown. This figure demonstrates that as a 
result of the stripes rj =^ for all nuclei on the boundary 
between stripes, despite the fact that the underlying lat- 
tice is square. It should be strongly emphasized, however, 
that we use a stripes-based toy model for the demonstra- 
tion because of its simplicity, and our experiment has no 
bearing on the model developed in Ref. lQj where the 
stripes are dynamic on a time scale much shorter than 
our experimental one. 

Our objective in this work is to determine 77 as a func- 
tion of doping in different environments using ADNQR, 
which is described in detail elsewhere |ll| . In this modifi- 
cation of standard NQR one measures the echo intensity 
in a coil, whose symmetry axis is in the direction 

f = (sin 9 cos 4>, sin 9 sin </>, cos 8) 



with respect to the principal axis of the EFG, as a func- 
tion of 9 and (p. This situation is demonstrated in the 
inset of Fig. The echo is obtained using & tt/2 — t — n 
pulse sequence. The pulses are of frequency /, duration 
t„/2 an d 2^/2 7 an d create a field Hi = W1/7, where 7 is 
the copper nuclei gyromagnetic ratio. The basic idea be- 
hind this experiment is that when r is in the z direction, 
and 77 = 0, no spin transitions will occur and no echo 
will be formed. In contrast when rj =/= an echo will be 
formed even though r is parallel to z. For an arbitrary 
9 and <p the nuclear magnetization M in the coil, at the 
time of the echo, is given by 

M(r),9,(f>) = M Asin 3 (Awit f ) (4) 

where 

A= — 1 {[Q + r^-t-firir-osm)] sin 2 (6>) + [2?7Cos(6l)] 2 } 1 
2V3 + r] 2 

(5) 

and Mq is a constant. 




e [degrees] 



FIG. 2: (Color online) The expected echo intensity for various 
values of r\ in an oriented powder (averaged azimuthal angle 
4>) as a function of the polar angle 9 between the RF field Hi 
and the z direction of the EFG. The pulse length is optimized 
at 9 — 90 and kept constant throughout the sample rotation. 
The inset shows the experimental configuration. 



We apply the ADNQR experiment to oriented powder 
of YBa2Cu30j, with c||z. In such powders the a and b 
directions are mixed. In this case M is obtained from 
Eq. 01 by averaging over cf>, namely, 

mm^^/m^m)^. (6) 
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This averaging must be done numerically, and the ex- 
pected M (77, 9) /M (77, 90) for various values of 77 is pre- 
sented in Fig [21 The pulse length (t n / 2 ) is optimized at 
9 = 90 and is kept constant throughout the rotation of 
the sample with respect to the coil. 

The ADNQR technique has three major strengths: (I) 
It allows the determination of 77 even for a wide NQR 
spectrum and at every point on the spectrum, thus pro- 
viding a great advantage over nuclear magnetic resonance 
(NMR), which can determine 77 only from the entire spec- 
trum with no local resolution; (II) unlike NMR, where 
one needs to fit the spectrum to 5 different parameters 
out of which only one is 77, ADNQR is sensitive only to 
77; (III) it allows the determination of rj without the ap- 
plication of a static magnetic field and can be used even 
in the superconducting state. The weakness of ADNQR 
is that it is very insensitive at small 77's as demonstrated 
in FiglU 
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FIG. 3: (Color online) NQR intensities versus frequency / 
normalized by f 2 in three different YBCO,, samples. The 
YBCC-6.678 sample is fully enriched with 63 Cu. YBCO7 and 
YBCO6.08 have natural abundance of 65 Cu and 63 Cu. Also 
shown are Gaussian fits indicating the contribution to the 
spectrum of different Cu(2) environments. The arrows show 
the frequencies where ADNQR is applied. 



Before presenting the ADNQR results it is essential to 
review the NQR frequency assignments in YBCO. In the 
(a), (b), and (c) panels of Fig. we show the NQR lines 
at a temperature of 100 K for three different samples with 
y = 7.0, 6.678, 6.08, respectively. The y = 6.678 sample 
is fully enriched with 63 Cu. The lines are obtained by 



a spectrometer with a home-made automated frequency 
sweep. The echo intensity I is normalized by / 2 , where / 
is the applied frequency, in order to correct for population 
difference and the induced signal in the coil. The solid 
lines in Fig.|21are fits of Gaussians to the data indicating 
that the spectrum is constructed from contributions of 
different local environments. 

It is simplest to review the enriched y — 6.678 spec- 
trum in panel (b) first. This spectrum exhibit three peaks 
at / = 28.1, 29.2, and 30.8 MHz (none of which are from 
the 65 Cu isotope). The peaks are classified in terms of the 
number of oxygen surrounding the chain copper [Cu(l)] 
neighboring the detected Cu(2). The higher the oxygen 
coordination of the this Cu(l), the higher the frequency. 
In other words, the peak in / = 30.8 MHz is from Cu(2) 
adjacent to a Cu(l) in a fully oxygenated environment 
[Cu(2) 4 ], the peak in / = 28.1 MHz is from Cu(2) ncaring 
a Cu(l) in a fully oxygen deficient environment [Cu(2)2], 
and finally, the peak in / = 29.2 MHz is from Cu(2) 
whose neighboring Cu(l) is missing one oxygen [Cu(2)s]. 

The main features in the spectrum of the natu- 
ral YBCO7 (without enrichment) depicted in panel 
(a) are two resonance peaks at / = 31.55 MHz and 
/ = 29.3 MHz. In addition, each peak has two 
shoulders to the left (31 MHz/28.7 MHz) and right 
(32 MHz/29.7 MHz) of the biggest peak. The origin of 
the right shoulder is not clear. The left shoulder is associ- 
ated again with a partially oxygen deficient environment 
(mostly Cu(2)a). This spectrum, including the shoul- 
ders, resembles the spectrum obtained by others 0,0], 
where the lines are identified as the plane-copper iso- 
topes [ 63 Cu(2)] and [ 65 Cu(2)]. The ratio between the 
peak intensity of the two isotopes is the one expected 
from their natural abundance, indicating that our spec- 
trometer functions properly. 

It is important to mention that the Cu(l) can also con- 
tribute intensity to the Cu(2) signal 0|. This is demon- 
strated in panel (c) for YBCOg.os- In this antiferromag- 
netic sample the Cu(2) signal is at / = 90 MHz 0| and 
in panel (c) only the two isotopes of Cu(l) can be seen 
|l4j . This panel clarifies our careful choice of oxygenation 
level of the underdoped compound. The Cu(2) peaks of 
the 77 = 6.678 sample are not contaminated by the Cu(l) 
signal. 

The ADNQR experiments in the YBCO7 sample were 
done on the 63 Cu(2) peak and the two shoulders. In 
YBC06.678 we investigated the peaks of all different lo- 
cal environments. The frequencies where ADNQR was 
applied are marked by arrows in panels (a) and (b). The 
ADNQR experiments were also done at T = 100 K tem- 
perature using an automated sample rotor. To improve 
Hi homogeneity we used a spherical coil. 

The ADNQR results for YBCO7 and YBC0 6 . 678 in the 
9 = to 180 range are depicted in Fig. 0] In all cases the 
intensities at 9 = and 180 are lower than at 9 = 90, as 
it should be. It is also clear that the intensity as a func- 
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tion of angular deviation from 8 = 90 drops faster for the 
Cu(2)4 in both dopings. This is a model-independent ob- 
servation of the fact that rj is smallest for Cu(2) 4 . The 
fit of Eq. Elto the experimental data is demonstrated by 
the solid lines. In the fit we allow a common finite base 
line for a given sample, in order to account for some un- 
known amount of misalignment. We found 24% and 17% 
misaligned crystallines in the YBCO7 and YBC06.678j 
respectively. This misalignment is typical. The best fit 
for YBCO7 peak gives ij — ± 0.15 in agreement with 
values previously obtained by NMR for optimally doped 
YBCO, where the values r\ = and rj = 0.01 ±0.01 
were found. This result reassures us that the experimen- 
tal method and analysis are correct. Surprisingly, our 
fit also gives 77 = ± 0.15 for the Cu(2) 4 in YBC0 6 . 6 78 
sample (/ = 30.8 MHz). In contrast, in oxygen deficient 
environments in both YBCO7 and YBCOg.678 V — 0.6. 
This is the major finding of our work. 
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FIG. 4: (Color online) The echo intensity as a function of the 
angle 9 between the RF field Hi and the principal axis of 
the electric field gradient z for the two YBCO samples with 
different oxygen level y, and at three different points in the 
specturm shown in Fig. EH a) and (b). The solid lines are fits 
to Eq. E]as described in the text. 



The data from YBCOg.678 are the first evidence for 
charge anisotropy in such highly doped YBCO. As men- 
tioned before, to date, static charge/spin inhomogeneity 
has been detected in YBCO by neutron scattering and 
//SR only up to y < 6.35 However, the possibility 

of electronic phase separation, reflected by 77 > 0, is found 



(within our experimental sensitivity) only in conjunction 
with quenched disorder. This means that Cu(2) 4 cannot 
reside on a static boundary between two phases of dif- 
ferent charge concentration, like the boundary generated 
by our toy model. This result is in accordance with re- 
cent STM measurements 0. We thus conclude that in 
YBCO charge inhomogeneity stems from quenched dis- 
order. 
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